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I. INTRODUCTION 

A great deal of attention has been focused on the studies of the electronic structure of 
the high T^ superconducting cuprates (HTSC) in the past decade [|l]-|^. In the undoped 
case, the charge mobility is suppressed by the strong on-site Coulomb repulsion and the 
spin-i antiferromagnetic Heisenberg model is believed appropriate to describe the electronic 
properties of these materials [Q. In the opposite limit, when the cuprates are overdoped in 
the copper-oxide plane, the HTSC is metallic and the observed electronic structure is in good 
agreement with predictions of the local density approximation (LDA) calculations ||^ . In- 
between these two limits, many interesting phenomena appear including superconductivity. 
The evolution of the electronic structure from one limit to another should play a key role 
for the understanding of numerous anomalous properties of HTSC discovered in this doping 
range. 

When slightly doped with holes in the antiferromagnetic Cu-0 plane, the Fermi surface of 
the mobile holes is expected to be small pockets centered at (±7r/2, ±7r/2). On the contrary, 
in the overdoped case the Fermi surface should be a large one centered at (tt, tt) as reported 
by both angle-resolved photo-emission spectroscopy (ARPES) experiments and the LDA 
calculations [Q. However, the crossover from one limit to another is not well understood 
yet. 

The recent ARPES experiments [H-§] on underdoped samples seem to indicate that 
the Fermi surface of these compounds is probably half-pocket-like, i.e. a small pocket near 
E (7r/2, tt/2), but lacking its outer part in the reduced Brillioun zone scheme ( see the dashed 
line in Fig. 3). The data of Marshall et al. PJ^ show a well defined Fermi- level crossing 
along the (0, 0) to (tt, tt) direction, but no such crossing has been detected along the (vr, 
0) to (tt, tt) direction. This result differs substantially from that of the optimally doped 
samples which exhibit a large Fermi surface with well defined Fermi-level crossing in both 
(0, 0) to (tt, tt) and (vr, 0) to (vr, tt) directions. 

Another puzzling phenomenon in underdoped cuprates is the presence of a pseudo-gap 



in the normal state. Both Stanford |^ and Argonne [^ groups have shown strong evidence 
for the suppression of the density of states near the Fermi level indicating at the formation of 
the pseudo-gap. Several different measurements, including specific heat and thermal power 
PJTO[], Hall coefficient [jll|, and transport ||T2[, as well as nuclear magnetic resonance (NMR) 



|T^ and neutron scattering experiment |T^ , have also shown the presence of the pseudo-gap 
in underdoped cuprates. The symmetry of the pseudo-gap has been found d-wave like |0,§] 
with a maximum value in the (vr, 0) to (tt, tt) direction, and it is very much reduced in the 
(0, 0) to (vr, vr) direction. This result is also consistent with earlier data of Marshall et al. 
0, where dramatical changes have been observed for the spectral density near (0, tt). 

Theoretically, the origin of the pseudo-gap can be ascribed either to short range AF 
fluctuations in the weak coupling approach ||15| or to d-wave pairing of spinous ||T6| , p!7| in 



the resonant valence bond (RVB) state [T^. Due to the SU(2) symmetry at half filling, the 



d-wave and s-l-id wave [|19| RVB states, as well as vr flux state pO| are degenerate ||2T| . Away 



from half filling, the commensurate flux phase (CFP) |^, the staggered flux phase (SFP) 
PB[ and the d-wave RVB state [^ are competing candidates for the ground state. Further 



study shows that the CFS is favoured only in the limit t <^ J [^, which is not the actual 
case in cuprates. Zhang P3 has studied the competition and possible coexistence of the SFP 



and d-wave RVB state, and he concluded that SFP is unstable against d-wave RVB state 
and there is no coexistence of these two states at any non-half-flllings. However, Zhang's 
conclusion was based on the renormalized mean-fleld theory which satisfies the no-double 
occupancy constraint only on average. To improve the situation, in one of our earlier papers 



25| , D.N. Sheng and two of us (Z.B. Su and L. Yu) have reconsidered this issue using a 
better treatment of the slave-boson mean field theory, proposed by Nori, Abrahams and 
Zimanyi (NAZ) |2^, transmuting the holons into spinless Fermions via a 2D Jordan- Wigner 
transformation. Using NAZ statistics transmutation trick it was found that the vr-fiux state 
can coexist with the s-|-id-wave RVB state below a critical doping 6c, while the d-wave state 
takes over beyond that critical concentration. 



Recently, a SU(2) slave boson theory was proposed by Wen and Lee [jl3, significantly 



improving the earlier f/(l) version, especially for the underdoped systems. In their approach 
the SU(2) symmetry is preserved away from half-filling by introducing two kinds of slave 
bosons. Considering the effective interaction between holons and spinons they obtain a 
small Fermi Surface in underdoped regime. We will elaborate further on their results later, 
comparing them with ours. 

In the present paper we develop a slave boson approach for the t — J model, with a 
more careful treatment of the hard-core constraint for the holons by transforming them into 
fermions. Following NAZ |^, a statistical gauge field Aij is introduced to implement this 



transmutation. The effective gauge field seen by the holons will be a superposition of the 
TT-fiux contributed by the spinons and the statistical gauge field reflecting the hard-core 
nature of the holons. Using the equivalence of the vr-staggered and vr-uniform flux, we find 
that the statistical gauge field can substantially compensate the 7r-fiux. Then the effective 
flux acting on holons satisfies the commensurability condition which minimizes the holon 
free energy. The mean field phase diagram is obtained and the coexisting d-wave RVB and 
TT-fiux state is found energetically more favourable with respect to the pure d-wave sate in 
the underdoped region. This confirms the earlier result [^, and disagrees with Zhang's 



conclusion |Q, excluding the coexistence of these two states. It appears to us that the 
main reason of this discrepancy is due to the fact that in Zhang's approach the holons and 
spinons are not treated on an equal footing. The dynamics of holons is basically ignored in 
his approach, leading only to the renormalization of the hopping integral. On the contrary, 
in our approach, the full dynamics of holons is considered and the holons and spinons are 
treated equally. ( We will elaborate more on the differences of these two approaches in 
the concluding remarks.) Furthermore, we have calculated the physical electron spectral 
function in the low doping range and have found the half-pocket-like Fermi surface as well 
as the pseudo-gap in the same doping range. We emphasize that the pseudo-gap obtained in 
our approach opens right at the Fermi surface and therefore it is consistent with the ARPES 
experiments. Thus our calculation can very well explain the APRES results without invoking 
the additional (not physical) SU(2) symmetry away from half-filling. 



The rest of the paper is organized as follows: The mean field theory is given in Sec. II, 
whereas the spectral function for the physical electron is presented in Sec. III. Finally, we 
make a few concluding remarks in Sec. IV. 



II. MEAN FIELD THEORY OF THE t - J MODEL 

We begin with the slave boson representation of the t — J model in which the physical 
electron destruction operator Qo- is represented in terms of the pseudo-fermion operator fia 
and the slave boson operator hf as Qo- = fiabt- The partition function reads as : 

Z = JllD[f^Ar)]D[fUr)]D[bt{r)]D[k{r)] exp[- f Ldr] (1) 

with 

L = J2f,Ur){-^-f^2)Ur) + T.K{r){-^-f^iMr)+Nf,,6+N^,2{l-S) 

ia i 

-\ E(^^. + i bfb,) (r,. + 1 bp,) - ^ 5: btkb^b, + ^ 5: bfb, 

{ij) (ij) « 

-i5:A+.A,,+t'5:/+/,v6A^, (2) 

{ij) (ii') 

where we define 

(T (7 

The chemical potentials for the boson fii and the fermion /i2 are chosen to satisfy the global 
constraints Uf + Ub =1 and nb = 6, with 6 as the doping concentration. We take the 
exchange integral J = 1, while t, t' are the nearest neighbor, and next nearest hopping 
terms, respectively. Hereafter (ij) means nearest neighbor, while (ii') means next nearest 
summation. A^ is the total number of sites, and f3 is the inverse temperature. 
We introduce the Hubbard-Stratonovich (HS) transformation: 
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^= / ]lD[ft.]D[UD[bt]D[hmD[xUD[x^mD:mD^. 



{ij) 



■\{D[pl,]D[pu^] eM' fL'{r)dr] 



(3) 



with 



ia i {ij) 



(ij) (ii') {ij)a 

-In [Dh{r)A.j + h.c] - I Y.(xl,hth, + h.c.) + 1' ^ [p*,/+/.. + h.c. 



{ij)<^ 



{ii') 



{ii')(T 



-^^izbthbp, + ^^J:Kb^+ piT.^ + p^Ei^ - ^)- 



(ij) 



2 ^ 



(4) 



In the saddle point approximation of the above action, we treat the amphtudes of the 
HS fields D^, Xij and pu/ as site- and time-independent with some phase factors, namely, 

Dijir) = A r]ij, Xiji^) = p exp{i(j)ij), pu^ir) = p 



Vij 



with 

1 (j = i±x), 

-1 (j = i±y)., 

while (pij is the phase factor in the Affleck-Marston flux phase 
later. 

We then obtain the mean-field Hamiltonian per site as: 



to be described in detail 



1 



1 



Hmf = Hc + Hs + Hfi, 



t' 



H, = ^A' + -/ + /X2(l -5)+ p,5 - -5(5 - 1) + t'\p\\ 



A 



Hs = -^i: (e^'^-7+/,. + h.c.) - ^ E h^/^/;. + h.c 



2N 



foV 



2N 



feV 



+^ H (/ii/iv + h.c.) -^Yl ftfi^ 

(ii')(T icj 

Hh = -^ Ue'^'^Kb, + h.c.) + ^ El^.'-^^' + ^-c-) - ^ E &.^&^- (5) 

We introduce a statistical gauge field operator Aij to account for the hard-core nature 
of the holons, namely, bfbj = exp{iAij)hfhj, with hf, hj as fermion operators. In the mean 
field approximation, the statistical field A^j is considered as a c-number, but it must obey 
the constraint: 

E Aij = 27i{2q + 1)6 

plaquette 

with q as an integer.. This constraint can be satisfied by introducing an additional uniform 
flux attached to the holons. Then the effective mean field Hamiltonian for holons reads as : 



^P v^ rjf<Ai,+A-,U + L , L ,1 ^P 



2N 



/ii 



^^ = -^ E [e'^^^^^'^^^hth, + h.c] - ^ E \e'^^"'-''"'^hth, + h.c] - 9 E hfh.. (6) 



(ij) {ii') 

The effective flux that acts on the holons is then : 



N 



^f^= J2 {Aj + (pij) + SttA; = 27f6 + Airq6 + 2'Kk + vr, 

plaquette 

where k is an integer. \i 5 = {2pi — l)/4A;i, we can choose proper gi and ki to satisfy the com- 
mensurability condition $/i = 27r5 which minimizes the total free energy of holons . In the 

present calculation we choose ki = 25 and pi = 1, 2 , 25, so the above commensurability 

condition is always satisfied. Thus we transform the holon Hamiltonian into a simple form 
equivalent to the problem of noninteracting lattice fermions moving in a uniform magnetic 
field. We then diagonalize this mean field Hamiltonian of holons by solving the well-known 
Harper's equation. 

For the spinon Hamiltonian we choose a proper gauge shown in Fig. 1, where (pij = | 
if Rij is along the arrow shown in the figure and (pij = — ^ in the opposite case. Then we 
divide the original lattice into two sublattices, labeled by 1 and 2. The spinon Hamiltonian 
written in the /c-space is found to be: 

1 1 

^^ = -]^E'(e(^) + P2)flkaflka - ^E'(e(^) + P2)f2kaf2ka 
kcr ka 



- ^E'^ [(7+ + n-)fLf2k. + h.c] + ^'A7- [f^,,U-,, + h.c 
where Y^ means summation over the reduced Brillouin zone with reciprocal vector Q 



(7) 



(7r,7r), 



7± = (cos kx ± cos h 



y)i 



e(k) = 2t'p [cos(kx — ky) + cos{kx + ky)] 



The above Hamiltonian can be diagonahsed by the Bogohubov transformation, and we 

obtain the quasi-particle energy as ii^^s, and with s = ±1: 

^1/2 



E^ 



ks 



sp 



Icos (kx) + COS {ky 



H2 - e{k) 



2, ,2 



+ A^7. 



(8) 



Now we calculate the spinon and holon free energy for given values of p, 5, p, /ii and 
/i2. The parameters are chosen as t = 2, t' = 0.2, and the free energy is calculated at 
temperature T = 0.05. We minimize the total free energy by varying the order parameters 
p, A, p and fi to obtain the mean field phase diagram as a function of doping concentration. 
The spinon dispersion obtained by our approach shows that the pseudo-gap opens up right 
at the Fermi level which vanishes along the line k^ = ^ky reflecting the d-wave symmetry. 

As shown in Fig. 2, the coexistence of vr-flux and d-wave RVB state is favored in energy 
in the low doping range, while in the high doping range the pure d-wave RVB state takes 
over. The critical value of the doping concentration is found to be near 0.3. We believe this 
phase diagram is correct only qualitatively, since the fluctuations around the saddle point 
have been ignored in our present study. This calculation has confirmed the earlier result 



| 25| . The new ingredient here is to introduce the next nearest neighbor hopping t' which 
turns out to be essential for obtaining the correct shape of the spectral function. Like in 



the CFP [^, the mean field treatment of the commensurability condition $/j = 2t[5 for 
holons breaks the parity and time-reversal symmetry at arbitrary filling. However, this is 
an artifact due to the mean field treatment. 



III. ELECTRON SPECTRAL FUNCTION 



Now we calculate the spectral function for physical electrons. It is more convenient to 
use the boson representation for holons instead of the fermion representation because we 
can not fully restore the Bose statistical property of the holon operator in the mean field 
approximation when the operator Aij is treated as a c-number. The boson representation 
for holons is a further approximation compared with the mean field treatment in the pre- 
ceding Section. However, we do use the coexisting 7r-flux and d-wave RVB order parameters 
obtained there to calculate the electron spectral function. Presumably, the main physical 
consequences of the coexisting phase are still kept. 

First diagonalize the holon Hamiltonian in this boson representation. After Fourier 
transformation Eq.( |) reads as : 

^h = -T^Il'(e(A;) + f^i)btkbik - T^Il'(e(^) + /^i)&4&2fc - ir^Il'7^ [(7+ + h~)btkb2k + h.c. 

k k k ^ 

(9) 

where the indices 1,2 represent the two sublattices defined above. 
This Hamiltonian can be diagonalised as: 

^^ = ^E'EtiiMi^, - <«,) - ^^^^^E'(/5.^/5. + <«.) (10) 

q q 



with 

pt 

72 



Eb{,q) = —^\/cos'^{qx) + cos'^ (qy), 



which describes two kinds of free bosons labeled by a and (3. At low temperatures, bosons 
occupy the low energy states. So we need to consider only the bottom of the lower band 
(the a band) near F point (0, 0). Expanding Eb{q) near (0,0), we obtain the effective holon 
Hamiltonian as: 

Hh = Y.i^^-f^B)at(^q, (11) 

where fiB is the renormalized boson chemical potential, and 

1 tp 

-f + 2t'p 



2mB 8 
10 



The physical electron Green's function is: 

G{k,r) = -{TAcUr)cLm), (12) 



where 



Hence 



1 \/2 

Cka = —Y^Ciaexp{-ik ■ Ri) = —{Cika + C2ka)- 



G{k,T) = -W{T^{c^Ur)4kM)) 

= -^ E (rr(/..+,.(r)/;,+,.(0)))(r.(6+,(r)6^,(0)). (13) 

We see that the physical electron Green's function in the slave-boson approximation is 
nothing but the fermion-boson polarization function. The spectral weight is thus: 

Im G{k, uj) = -— Y, 'Sau{k + q)Sl^{k + q)Dl^{q)Dfii{q)5{uj + ujg- e^^^) [^^(^g) + ri/(e^+J 

afiqv 

(14) 
where S and D are the transformation matrices used by us to diagonalize Hg and H^, 



respectively. The first term of Eq. (|1^) contributes a quasiparticle like peak centered at 
e^—fiB, being called the coherent part. The second term represents a featureless incoherent 
part extended over a wide range. The density of occupied states probed by ARPES is 
linG{k,u!)nf{uj). If we consider only the holon states with low energies as discussed above, 
we find that the total spectral weight of the coherent part, detected by ARPES equals: 

— J2' J dujSc^^ik + q)S*p^{k + q)Dl^{q)Dpi{q)5{u + Uq- el^^)nB{u}q)nf{u) 

a/Bqu 

^ ^ E 'SaAk)S;MD:MDpmnBiuj,)nfiek) 

aPqv 

= ^T.'Sa.{k)S;MS Ufie,) = i < Jtfk. > 6 = ^nMS (15) 

Since the spectral weight itself is proportional to ns{k), we conclude that the shape of the 
Fermi surface is approximatively determined by the spinon momentum distribution ns{k). 
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The width of the peak is roughly proportional to T^/^, which is very similar to Lee and 
Nagaosa's [^ result for the uniform RVB phase. 

Fig. 4 presents the momentum distribution of spinous in different directions (as defined 
in Fig. 3) within the slave boson approach. We see that the curve a which is along the 
direction (1,1), i.e. from (0, 0) to (vr,7r), shows a clear step in the momentum distribution 
indicating at the Fermi-level crossing. However, away from the (1,1) direction, as shown by 
curves 71 and 72, the step in the momentum distribution is smeared out by d-wave RVB 
order parameter, and the shadow band (the second peak beyond the Fermi surface) appears. 
Curve /5 shows the momentum distribution along the direction (7r,0) to (7r,7r), and there 
is no step detected, indicating the absence of Fermi-level crossing in that direction. The 
pseudo Fermi surface is plotted in Fig. 3 for 5=0.2, appearing as an unclosed small pocket. 

We then calculate the physical electron spectral function in the a, /3, 71, 72 four directions. 
The typical total spectral weight of a and 71 cutes, including both coherent and incoherent 
parts, is shown in Figs. 5-6. Fig. 5 shows the spectral weight along the a cut ((1,1) 
direction). We see the entire coherent part moving across the Fermi level in this direction, 
and no signature of an energy gap (depletion of spectral weight at the Fermi level) is noticed. 
In Fig. 6, the pseudo-gap is visible at the Fermi level, and the coherent part can still move 
entirely across the Fermi level. Along the 72 direction the pseudo-gap is bigger than what we 
find in Fig. 6., which reflects the d-wave symmetry. On the other hand, the spectral weight 
does not drop to zero abruptly, decreasing continuously as the k vector moves along the 72 
direction. Finally, for the /3 cut we can hardly see any coherent part due to the small value 
of the spinon distribution ns{k), and the spectral weight is pushed down to the range much 
lower than the Fermi level, indicating at the formation of a maximal gap in this direction. 
Again, no abrupt Fermi-level crossing is detected along this direction. Our result is fully 
consistent with the ARPES data obtained by Loeser et al., which show a nearly perfect 
Fermi-level crossing along the a cut and no Fermi level crossing along the (3 cut. They could 
also find the pseudo-gap from their experimental data, with vanishing value along the a cut 
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and a inaxiinuin along the /3 cut. Similar results were also reported by H. Ding et al. 0. 

In Wen and Lee's approach ||1^, the electron Fermi surface is determined by the pole 
of the renormalized electron Green's function due to interaction with the gauge field. An 
elliptical Fermi surface with strongly suppressed spectral weight on the outer edge is obtained 
in the low doping regime. This elliptic surface is not centered around S point, which is on 
the outer side and the wave function renormalization factor is exactly zero at that point. 
Hence there is no "shadow band" in their approach outside the boundary (0, vr) - (vr, 0). A 
spin gap opens outside the small elliptical Fermi surface, but the gap is identically zero at 
the entire Fermi surface. Also, the quasiparticle peak crosses the Fermi level in all directions 
spanned by the Fermi surface. The main features of our results are similar to theirs, namely, 
small half-pocket-like Fermi surface near E, and d-wave like gap in the excitation spectrum. 
However, there are several important differences. In our approach the spinon Fermi surface is 
adjusted by the chemical potential, i.e. the doping concentration, so the d-wave like pseudo- 
gap is "locked" on the Fermi surface, and the pseudo-gap vanishes only at four points at 
which the Fermi surface intersects the line k^ = ±ky. As soon as the points of the Fermi 
surface slightly deviate from the (1,±1) direction, a finite energy gap appears. Also, the 
quasiparticle peak of the spectral function does not cross the Fermi level in any direction 
except for (1,±1). Another difference is we have small pockets centered around S and we 
do have shadow bands beyond the line (0, vr) - (vr, 0). In principle, these different predictions 
can be checked by APRES experiments. It seems to us, however, that the presently available 
resolution is still not enough to resolve this issue. 

IV. CONCLUDING REMARKS 

Using the slave-boson approach, within the mean field approximation, we have shown 
in this paper that coexisting vr-flux phase with d-wave RVB state is energetically more 
favorable in the low doping range (below 6c ), compared with the pure d-wave RVB state. 
We have also shown that this coexisting phase can provide a rather natural explanation for 
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the half-pocket-like small Fermi surface, a d-wave symmetry pseudo-gap, and very different 
behavior for the spectral function in the a and (5 directions, in underdoped HTSC The main 
reason for this outcome is we have been able to better treat the single-occupancy constraint 
by introducing the additional statistical gauge field. 

The advantage of a mean field approach compared with more sophisticated theories 
taking into account gauge field fluctuations is that it gives a more transparent understanding 
of the experimental results. When studying the problem of a single hole moving in the 
background of antiferromagnetic fluctuations, an effective hole band is obtained by both 
numerical and analytic methods. The hole band minimum is located at (±|, ±|). Thus the 
Fermi surface of holes in the very low doping regime should be small pockets centered at 
(±|, ±|) due to antiferromagnetic fluctuations. This conclusion is confirmed by the recent 
ARPES results. Hence the antiferromagnetic fluctuations leading to a bipartial lattice must 
be considered in the low doping regime. In our approach we introduce the vr-flux phase 
which mimics the short range antiferromagnetic correlations which coexist with a weak d- 
wave RVB phase. Intuitively, the vr-flux phase carries the physics of half-pocket-like Fermi 
surface while the coexisting d-wave phase is responsible for the pseudo-gap opened at the 
Fermi surface. 



Zhang has concluded pj], however, that the vr-flux state is unstable against the pure 
d-wave RVB state, excluding thus the possibility of their coexistence, using a renormalized 
mean field theory within the slave boson approach. What is the origin of the discrepancy be- 
tween our result and his conclusion? We find that in the mean field theory of the t — J model 
the total energy is divided into two parts. One is the kinetic energy which is characterized 
by the normal average terms < /^/ > (for spinon) and < h^h > (for holons). The other 
is the exchange energy which is characterized by the anomalous average term < f^f^ >. 
The vr-flux state favors the kinetic energy, while the d-wave RVB state favors the exchange 
energy. In Zhang's approach the dynamical and statistical properties of holons which are 
shown to be nontrivial by our present study, are not considered, and the contribution of 
< h^h > is ignored. This underestimates the kinetic energy gained by the delocalized mo- 

14 



tion of holons. Therefore the coexistence of vr flux and the d-wave RVB state is excluded 
in Zhang's calculation. In our present study the dynamics of holons is fully considered, so 
we find that in the low doping regime the tt-Aux state coexisting with d-wave phase is more 
stable against both the pure 7r-fiux state and the pure d-wave RVB state. The admixture of 
a small amount of the d-wave state will provide holons with more energy gain at the Fermi 
surface compared to the pure vr-flux phase. Due to the smallness of the Fermi surface in 
the coexisting state the d-wave order parameter (of the order 0.1) is much smaller than the 
vr-fiux order parameter (of the order 1). The energy gap opens at the Fermi level due to the 
presence of the d-wave RVB order parameter. Since there is no Bose condensation in our 
calculation, this mixed state is responsible for the pseudo-gap in the normal state. With 
the increasing of the doping concentration, the area enclosed by the spinon Fermi surface is 
enlarged and the states near the Fermi surface take a higher percentage of the total occupied 
states and play a more important role compared with the low doping case. Therefore, in the 
high doping range, the d-wave RVB state gains free energy more effectively, and a crossover 
from the coexistence state to a pure d-wave state takes place. We believe our mean field 
treatment can serve as a good starting point for further improvement. It turns out that the 
mean field approximation of the U{1) x SU{2) Chern-Simons bosonization approach for the 
2D t — J model |^0|] leads to very similar results as the present study. The condensation of 



holons will lead to superconductivity which was not considered in this work. The fluctua- 
tions around the saddle point will create an effective attractive interaction between spinous 



and holons, as recently discussed by Wen and Lee |jT^ • The role of these fluctuations in our 
approach require further studies. 

Very recently. Ding et al. p8[ reported that there seems to be a large Fermi surface for 
HTSC over the entire range of doping, from underdoped to overdoped cases. However, their 
resolution is still not satisfactory for the substantially underdoped samples {Tc = 15K ). It is 
yet to be seen how strong is their evidence. Also, the transport and optical data have shown 



rather strong evidence in favor of small hole pockets in the highly underdoped samples PU 
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Figure Captions 

Figure 1. Diagrammatic illustration of the tt-Aux state. A 7r(— tt) magnetic flux flows 
alternatively through neighboring plaquette of the lattice. 

Figure 2. The comparison of the total free energy F for the coexisting 7r-flux + d-wave 
RVB state and the pure d-wave RVB state. 

Figure 3. The unclosed pocket-like spinon Fermi surface at 5=0.2. The four directions 
along which we calculate the physical electron spectral function are marked in the figure as 

a, /5,7i and 72. 

Figure 4. The momentum distribution of the spinous along the four directions defined in 
Fig.3. 

Figure 5. The physical electron spectral function along the a cut, from which we see a 
clear Fermi-level crossing and no pseudo-gap. 

Figure 6. The physical electron spectral function along the 71 cut. A pseudo-gap opens, 
but the Fermi-level crossing is still visible in this direction. 
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